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Abstract 



o 

u. 

f-n . It is observed that one of Einstein-Friedmann's equations has formally the aspect 

Qh| of a Sturm-Liouville problem, and that the cosmological constant, A, plays thereby the 

role of spectral parameter (what hints to its connection with the Casimir effect). The 
subsequent formulation of appropriate boundary conditions leads to a set of admissible 
values for A, considered as eigenvalues of the corresponding linear operator. Simplest 
boundary conditions are assumed, namely that the eigenfunctions belong to L 2 space, 
with the result that, when all energy conditions are satisfied, they yield a discrete spec- 
trum for A > and a continuous one for A < 0. A very interesting situation is seen to 
occur when the discrete spectrum contains only one point: then, there is the possibility 
■ to obtain appropriate cosmological conditions without invoking the anthropic princi- 

ple. This possibility is shown to be realized in cyclic cosmological models, provided the 
£Sj \ potential of the matter field is similar to the potential of the scalar field. The dynamics 

of the universe in this case contains a sudden future singularity. 

(N 



PACS: 04.20.-q, 98.80.-k 



^ ; 1 Introduction 

In Einsteinian Gravity, when quantum corrections are taken into account the determination 
of the cosmological constant (cc), A, is very much related with vacuum fluctuations and, 
eventually, the Casimir effect [1], in one way or other [2]. The problem of the cc is one 
of the most important and intriguing in modern cosmology and QFT. The traditional way 
to understand why A is so small (A < 1CU 120 in Planck units) is ordinarily based on the 
anthropic principle j3]. Up to now one has no reasonable alternative way to explain this fact 
in Einsteinian Gravity, unless one moves on to extended models with dynamical dark energy 
or modified gravities [I]. 
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Recently Barrow and Shaw j5] suggested a non-anthropic solution to the problem: by 
making A into a field and restricting the variations of the action with respect to it by 
causality, they managed to obtain an additional Einstein constraint equation. One can say 
that this approach is based on a different interpretation of the cosmological constant, notably 
to consider the cosmological constant to be a field variable. 

In the present work we use a somehow similar procedure to calculate the value of A, 
namely to interpret it as an eigenvalue of a Sturm-Liouville problem, rather than as an 
integration constant. Our approach is in fact more conservative than the method in [5] 
since we use the standard Einstein-Friedmann equations without any additional constraint 
equation. Instead, we consider one of Friedmann's equations as a spectral problem and look 
for a class of boundary conditions which may allow us to actually calculate the corresponding 
eigenvalues (in this way, our procedure is closely related to standard investigations of the 
Casimir effect [El d E] by spectral methods). A most simple condition is to impose that 
eigenfunctions be elements of L 2 . As we will show, this choice results in the following 
interesting consequence: if the universe is filled up with a matter field, <ft (except for A), 
such that all the strong energy conditions are satisfied, then one gets a point spectrum for 
positive values of the cosmological constant, and a continuous spectrum for negative ones. 

The next step will be to consider models which actually result in a positive discrete spec- 
trum with just one single point. This idea is inspired in the conceptual approach developed 
in the paper by Linde and Vanchurin |7j: one can fix all the parameters in the landscape, 
including the value of the cosmological constant, without using any anthropic considerations, 
what "will give us a chance to return to Einstein's dream of a final theory, which may allow 
us to make sharp and unambiguous predictions despite the abundance of choices". Clearly, 
this case of a one-point spectrum is a good example of a possible realization of this cherished 
dream. 

An additional bonus of our approach is the surprising connection that will be uncovered 
between the Steinhardt-Turok cyclic models |8J and Barrow's sudden singularity [9j— or sin- 
gularity of type II according to the classification in [10]. In short, on one hand, models with 
Steinhardt-Turok-like cyclic potentials may indeed have much to do with those singularities, 
and on the other hand, such models exhibit a single-point positive spectrum for A in the 
framework of our approach. This striking connection will be a crucial issue in this article. 



2 Friedmann's equation as a Sturm-Liouville problem 

A method for constructing and analyzing exact cosmological solutions of the Einstein equa- 
tions based on representing them as a second-order linear equation (we call this the lineariza- 
tion method, in what follows) was discussed in [TT] (other methods for constructing exact 
solutions in cosmology can be found in [T2]-[TT]) . Indeed, it is easy to see that, in the case of 
the flat Friedmann metric, the third power of the scale factor ip = a 3 satisfies the equation 

where p is the density and p the pressure of the matter filling the universe. In this paper 
we use the system of units with 8nG/3 = c = 1. In the case when a minimally coupled 
scalar field, 0, with the self-interaction potential V(<f)), is dominant, and in the presence 
of a cosmological constant with density A, Eq. ([I]) formally coincides with the Schrodinger 
equation 

§ = (U-W, (2) 
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where the potential is U(t) = 9V, and the spectral parameter is A = — 9A. In (j2J), the 
quantity V is assumed to be a function of time: V(t) = V(4>(t)), which has been called the 
history of the potential in [UJ. Giving an explicit form U (t) together with the corresponding 
boundary conditions allows to find the general solution of Eq. (J2J). An important conse- 
quence of this investigation is that the regime is independent of, or just weakly dependent 
on the type of the potential, what is quite significant for the whole theory. Unfortunately, 
the ubiquitous problem of the end of inflation turns out to be substantially more difficult 
here, and solving it in the framework of the approach described above apparently involves 
additional assumptions. (The authors of proposed modifying the potentials to make 
them depend on the temperature. Stablishing a Friedmann regime can then be described 
as a phase transition in the matter state of the early universe) . The study of Eq. in its 
applications to cosmology was continued in |18j and |19| . where the Darboux transformation 
was used to construct new exact solutions (a similar technique was used in [20] to obtain 
exact solutions on the brane and on the encompassing space carrying an orbifold structure). 
We consider the Einstein equations in the Friedmann metric, 



a 2 



a* 



and assume that a = a(t), p 



a 
a 

P(t), P 



k 

P ~ — , 
a z 

-\{P + Zp) 



(3) 
(4) 



p(t) is a solution of these equations, for k = 0. 



Then, the function ip Tl 



where the "potential" is 



is a solution of the Schrodinger equation 



= U n (t)lp r , 



(5) 



U n {t) 



n 2 p 



3tz 

y(p + P)- 



(6) 



If the universe is filled up with a scalar field with the Lagrangian L 

n(n — 3) 



^■-V (</>), then 



-<j) 2 +n 2 V((j)). 



(7) 



The effectiveness of the method presented in [llj (and developed in [IB] and |19|). pre- 
cisely consists in reducing a complex nonlinear problem into a linear equation. This allows to 
obtain a complete set of two-parameter solutions, which exhibit inflationary behavior under 
very general assumptions. The fact that U basically coincides with V was not used anywhere 
in those papers and consequently played no role there. In a similar way, we here consider a 
generalization of this method to arbitrary n. 

Further to this point, the physical meaning of the "potential" U3 is clear only for a 
universe filled with a scalar field. Indeed, if we consider a universe in which, for instance, 
electromagnetic radiation is dominant, the physical meaning of the quantity C/3 becomes 
again ambiguous. 

We note that if a solution of (J2J) is known, then we can use Eqs. ([3]) and (J2J) to find the 
scalar field 



and the potential 



<j>{t) 



V(t) ^ 3 



V2 



dU 




(8) 



(9) 
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In principle we can obtain the dependence V = V((p) from these expressions, although it is 
clearly not always possible to do this explicitly. For example, if n = 1 and 



U(t) = At 2 [At - 3 
then 

p=16t 6 , p = 8t 2 (l-2t 4 ), w = -l + _L H = -4t 3 , 

and 

V($) =4V2$ 3 + 4 v / 3$ 2 , 

after performing the Bogoliubov transformation <p = l/\/6 + $. 
In the general case, a solution of fl5]) has the form 

*n = Cl-0n + C 2 ^n, (10) 

where ^„ is a linearly independent solution with the same potential: 

rt df 

11) 



Equation ( [101) allows to prove the following statement. 

Assertion. Let a = a(t) be a solution of (J3J) and (jl]) for /c = and the corresponding p and 
p. Then, the three-parameter function a n = a(t; Cx,C2,n), of the form 

( f dt \ Vn 

a n = a [ ci + c 2 J — J , (12) 

is a solution of Eqs. (|3]) and (|1J) for the new energy density p n and pressure p n satisfying the 
condition 

n Pn - Y (Pn +Pn) = U p ~ — {p + p) . (13) 

Finally, if we assume that the universe contains a nonzero vacuum energy with density 
p A c 2 , in addition to the matter fields, then Eq. (j3J) takes the form of the spectral problem 

= (U n (t) - A n )^„, (14) 

where the spectral parameter is A n = —n 2 p A . Just as for Eq. (T5]) , we can consider a problem 
for the eigenvalues and eigenfunctions of Eq. (I14p . if we specify homogeneous initial condi- 
tions. As noted in [TT], Eq. (j2j) has the form of a quantum mechanical problem with a discrete 
spectrum. The fact that each solution of this kind only admits a bounded or countable set 
of allowed values of the cosmological constant (if we specify homogeneous initial conditions) 
may help a lot to clarify the question of the actual value of the cosmological constant. It is 
clear that, in the case n = 1, the SEC results in the condition U(t) < 0. If one would like to 
obtain the discrete spectrum of A for such potentials, one would need to suggest U(t) — > 
at \t\ — > oo and this readily means that one can get the discrete spectrum for positive values 
of A only. 

This leads one to ask the following question: Is it possible to obtain models with a 
discrete spectrum containing just one single point? 

For answer one need to use some boundary condition on ip n . In this article we suggested 
that ip n is element of square-integrable functional space L 2 . This is just the hypotheses 
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about boundary condition. Therefore, our aim looks as follows: to investigate the class of 
cosmologies which results in just one allowed value of vacuum energy, on conditions that 
ip n is element of L 2 . We'd like to understand rather physical ground which may results in 
this picture. In the issue, we have obtained two results, two cosmologies which admits the 
possibility under consideration: 

1) The universe with some new form of sudden future singularity (Sec. 3). 

2) The universe without exact sudden future singularity (rather with "smoothed singu- 
larity") but filled with scalar field <p with deep narrow self-acting potential V(<fr) (Sec. 4). 
Such potential is very similar with cyclic potential of Steinhardt-Turok (ST) except that the 
ST cyclic potential V(<p) has a different asymptotic behavior at — > ±oo. Therefore, to 
obtain ST potential (and breathe the life into our formal model!) one need to modified the 
model from the Sec. 4. We do it in the Sec. 5. 

Thus, we have two cosmologies which allows one to obtain the extreme form of non- 
anthropic solution of the problem of cosmo logical constant. Remarkably, these cosmologies 
are not absolutely new, speculative models; vice versa, these cosmologies are connected with 
two well known and popular (at present) models; Barrow sudden future singularities and 
Steinhardt-Turok cyclic models. And last, but not least: these two models are connected to 
each other by the following way - the sudden model (Sec. 3) may be obtained from the "ST" 
models by the limiting process. 



3 New class of sudden future singularities 

We assume that the matter field contains baryons, that there is also radiation and dark 
matter, and that the dark matter leads to a sudden future singularity at t — t s , namely 

p s (t s ) = p s = +oo, p s (t s ) = p s < oo, < a(t s ) < oo. (15) 

For the realization (TTS]) . let us choose the dark matter pressure to be of the form: 

2 

PDM{t) = Pdm (t) + -a S(t s - t). (16) 

Hence, pDAi{t) is the "regular" part of the function pDM{t), while a 2 is a positive constant. 
From the equation of state for the dark matter, pdm = wdmPdm, in which wdm = 0, t ^ t s , 
it follows that pDMif) = 0. Then, for the energy density and pressure, we have 

2 

p = \ 2 + p DM , p= -A 2 + -a 2 5(t s -t). (17) 

The contribution of the dark matter decreases with time for t < t s . For t = t , Pdm/^ 2 < 0.5 
and we can neglect Pdm in the equations. Therefore, Eq. (fT4^) can be written in the form 

(A 2 -a 2 5(t s -t))^ (18) 

Let us now assume that the function ip is continuous at the point t = t s . The first derivative 
of the scale factor has a gap at the time of occurrence of the sudden future singularity, and 



,(+) _ hi-) 



-a 2 a s , (19) 



where = ^ lim^d(t) and a s = a(t s ). From the physical point of view, this gap corresponds 
to a density jump in the singularity. One gets, 

Sp = di+)2 - = a2 ^ +) +^) , (20) 
a 2 a s 
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For simplicity one can set Sp = 0. Therefore, if the universe expands before the time t = t s , 
then at time t = t s contraction begins. It is clear that the solutions for t < t s and t > t s 
cannot be connected. But it is interesting to investigate this model further. 
The condition f )19p at 5p = leads to 



H 



s 
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Therefore, the density at the time of SFS is given by p s = a 4 /4 + k/a 2 s . One notes that for 
the open universe the following condition should be satisfied 

a 2 - 2a; 1 > 0. (21) 

Solving (I18|) with the conditions 5a s = and 5p s = yields 

a 2 a 2 
a-(t) = a s exp — (t - t s ), a + (t) = a s exp — (t s - t). (22) 

Therefore, the value of A is unique and, for the vacuum energy density, we have 

A=^,. (23) 

Condition (12T|) can be written as follows 

9 1 

s A 

4 Smoothed SFS 

Strictly speaking, the model in the previous Section has no physical sense. This is because one 
had to consider two different universes (with positive and negative Hubble roots, respectively) 
separated by the singularity. One needs to involve a smoothed "singularity", such that 
\p\ ~^ \Ps\ ^ 1; but \p s \ < oo and p s < oo. Let us take 

u -® ^ -t^fcry (24) 



If K — > oo, we have 



and taking into account that 



one concludes that 



/+oo 
U K (t)dt = -a' 



lim U K (t) = -a 2 8{t). 



The potential as a function of time, for various values of k, is depicted on Fig. 1. Therefore, 
for k, = oo the potential (124")) leads to a sudden future singularity, as considered above. We 
assume that k » but k < +oo. The solution of (TI~4"]) . for the potential ( 1241) . can be 
written in the form 

m = (i - er /2 F L - s , e + 8 + 1, c + 1, ^] . (26) 
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Here F is the hypergeometric function, and we have used the following definitions: 



, / n A 1 / , .2d; 2 
4 = tanh(Kt), e = -, s = - — 1 + y H 

The condition k — )■ ±oo corresponds to £ — )■ ±1. Finiteness of the function requires that 

e-s= -N, N = 0, 1, . . . 

and, therefore, for the spectrum 



2a 2 



A « = 2 lV 1 + ~ ~ ' 1 ~ h ' X - (27) 

For k — y oo, one can easily show that 

\ N = a 2 /2 -kN^ -oo, if N ^0. 

But there is the condition that e > (as for e < the function f l26l) diverges at £ — > ±1) 
Therefore, A > and from Eq. (12711 one has 



1 + — - 1 I > 2N. (28) 



For sufficiently large k this inequality can be satisfied for N = only The remaining single 
eigenvalue is given by 



A* = | ( y 1 + ^- - 1 | , for k > a 2 /4, 

and, in this limit, one gets for the cosmological constant this same value. 
For the eigenfunctions 



m = (i-er /2 F (o^+M+i.i^i) = (i-e y" = cosher, £ = i u 



'l + ^-l 



(29) 



and the Hubble parameter is 



H = dln f {t) = A(«)tanh(/ct), A(«) = - h/l + — - 1 I . (30) 
at 2 \ V ft 



For the matter pressure and density, we get 



A(/t) 2 a 2 K + X 2 (K) 



ym cosh 2 (tfi;)' ^ m 3cosh 2 (/tt) 
The evolution of the scalar field can be obtained from the condition p m + p Tl 



l-2 



2 a 2 K-2X 2 U) . . .. 

>o ± — v arctan(exp(«;t)). (31) 

n V 3 
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t 



Figure 1: The scalar field potential as a function of time for various values of k. The depth 
of the well grows, while its width decreases, with k. 



For simplicity, we choose the "+" sign in (1311) and set 0o = 0. The potential of the scalar 
field is similar to the sine-Gordon potential: 



V^) = L m >(»), $ = -^ _-__ (32) 

The potential of scalar field as function of is depicted on Fig. 2. For large the scalar 
field changes sufficiently slowly. The corresponding potential at t = 0, for k » 1, has a 
deep well: 

y(t = 0) » -«V3. 



5 Asymmetric potential 

The potential considered in previous section is very similar with cyclic potential except that 
the Steinhardt-Turok cyclic potential V((fi) has a different asymptotic behavior at (j) — > ±oo. 
One can easily construct this potential starting from an appropriate U(t). Let us here 
investigate the case when the potential U(t) is taken with the form of a potential well 

U(t) =0, t< -t, 

U(t) = -U , -t < t < 0, 
U{t) =u , t> 0. 

Here Uq,Uq, and r are positive constants. One can perform the parametrization Uq = 
—A cosh 2 £, Mo = —A sinh 2 rj, where ^ and r\ are new constants, and write the solution ([Tip in 
the form 
1) t<-r: 

V>i = Ciexp( v /Z At); (33) 
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2) -r < t < 0: 

3) t > 0: 



■02 = C*2 sin(A/— A sinh£ ■ t + <5); 



(34) 
(35) 



V'a = C3 exp(— V — A cosh?] • t). 

Here Cj and 5 are constants. We consider the solutions which are finite at t — > ±00. 

The function ip and its first derivative should be continuous at the points t = t and t = 0. 
These conditions can be written as 



^ 2 (0) ^ 3 (0) 



M-r) M-t)' Mo) ^3(0)' 

The first condition yields tan(\/— A sinh £ • r — 6) = — sinh£, and the second tan<5 
Solving these equations for \J— A and 5, one gets 

sinh£(l + cosh 77) 



and 



tan(v— A sinh £ ■ r) 



tan 5 



sinh £ — cosh 77 
sinh£ 



sinh g 
cosh r] ' 



(36) 



(37) 



cosh 77 

We now take the limits £ — > 00, r — > 0, and assume that 

Uqt = —A cosh £ • r = er, 

where a = const. This limit corresponds to increasing the depth of the potential well at 
—t < t < and narrowing it so that its area remains constant. We also assume, for 
simplicity, that 77 << £. 

Then, the condition fl36l) can be written as 



tan z 



cosh 77 



-\z, 



a 



2aexp(-£)/ v /z A. 



9 



From the condition z « 1, that is tan ,2 ~ z, the corresponding one on the value of \J — A 
follows: 

v/^A = ^— . (38) 

1 + COSh 7] 

Thus, the value of the cosmological constant is unique and it is completely determined by 
the width and the depth of the well, and by the parameter i], too. 

Using Eqs. ([8]), (Q we can obtain the evolution of the scalar field <f)(t) and potential V((f>). 
For simplicity we assume that k = 0. For n = 3, from these equations we have, in the case 
of flat spacetime 

(i) t < -t: 

0! = 1O , V{(j>) = (39) 

(ii) -t < t < 0: 

-A sinh 2 £ 



-A sinh ft/2+5/2) |, V{4>) = (cosh(V6(0 2 - 2O )) - 2)+A, 

(40) 

(iii) t > 0: 

03 = 030, V{4>) = -Asinh 2 //. (41) 

In fact we have here used an unphysical sharp-cornered potential well. In a more realistic 
situation one needs use a potential well with "rounded edges". In this case, the potential of 
the scalar field V(0) becomes a smooth function of <fi. Now, let us consider Eqs. (I40j) and 




put that Uq >> then 5 = Pi/2. Since, at t — — r, <p = <p2o — y2?7o/3r/4 and, at t — 0, 
~ 020, on e must conclude that 

U(02O-50) = -y(l-t/or 2 ), y(0 2O ) = -^. 
Therefore, in the interval — r < t < the potential exhibits a very narrow and deep peak, of 



width 5(f) = ^2Uq/3t / '4 and depth —Uq/3. It is our conclusion that in a more realistic case 
one will obtain a potential for the scalar field similar to the potentials in cyclic models. 

It is interesting to consider the conditions (136|) and (137|) for a well with finite depth 
and width. The parameters £ and r] can be expressed through the spectral parameter, the 
depth of the well and the height of the barrier. The width of the well, r, can be written as 
t = (3/\/Uq, where /3 is a new parameter. Then, Eq. (136|) turns into 



I y/x - 1(1 + Jl + Uox/Uo) 

tantfyjl - l/x) = V ' \ (42) 

x - 1 - ^Uqx/Uq + 1 

where x = —Uq/X. For < < 7r/4, Eq. (|42|) has only one solution. The root exists 



provided the condition tan/3 > \Juq/Uq is fulfilled. 

For illustration, choose /3 = 0.5, Uq = 1, and uq = 0.2, in which case r = 0.5. The 
numerical solution ( l4*2l) yields x = 191.393 and, therefore, —A = 0.00522. The parameters £ 
and r] are equal to 3.319 and 2.522, respectively. From Eq. (136|) it follows that 5 = —1.146. 

6 Conclusion 

The main result of the present paper is the proof that, considering the cosmological constant 
A as an eigenvalue of a Sturm-Liouville problem allows one to obtain a set of admissible 
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values for A. In the special case of a situation with discrete spectrum containing only one 
point, there is the possibility to strictly determine the cosmological conditions leading to this 
value of A without recurring to the anthropic principle. We have shown in detail that this 
possibility is indeed realized in cyclic cosmological models when the potential of the matter 
field is similar to the potential of the scalar field. The dynamics of the universe in this case 
contains a sudden future singularity. As a byproduct we have also demonstrated that there 
is a remarkable connection between the Steinhardt-Turok cyclic models and sudden future 
singularity ones. 
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